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Runaway electron distributions are strongly anisotropic in velocity space. This 
anisotropy is a source of free energy that may destabilize electromagnetic waves 
through a resonant interaction between the waves and the energetic electrons. In this 
work we investigate the high-frequency electromagnetic waves that are destabilized 
by runaway electron beams when the electric field is close to the critical field for 
runaway acceleration. Using a runaway electron distribution appropriate for the 
near-critical case we calculate the linear instability growth rate of these waves and 
conclude that the obliquely propagating whistler waves are most unstable. We show 
that the frequencies, wave numbers and propagation angles of the most unstable 
waves depend strongly on the magnetic field. Taking into account collisional and 
convective damping of the waves, we determine the number density of runaways 
that is required to destabilize the waves and show its parametric dependences. 

I. INTRODUCTION 

Relativistic runaway electron populations have been frequently observed in various plas- 
mas, e.g. large tokamak disruptions electric discharges associated with thunderstorms 



2| and solar flares [3j. Runaway electrons are produced when the electric field is larger 
than a certain critical field (E c ), and the accelerating force overwhelms the friction for high 
energy electrons. The anisotropy of the runaway electron distribution can lead to destabi- 
lization of electromagnetic waves through wave-particle resonant interaction. Several studies 
have shown that the velocity anisotropy excites electromagnetic waves mainly through the 
anomalous Doppler resonance 0, [^J. Once the instability is triggered the distribution is 
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isotropized due to pitch-angle scattering. Previous work [4H6J has considered whistler wave 
instability driven by an anisotropic electron distribution as a possible cause for the observed 
magnetic field threshold for runaway generation in large tokamaks 0, S|- These calculations 
relied on a distribution function that was based on an approximate solution of the kinetic 
equation in the case when the electric field is well above the critical field, a>l, where 

E_ 4ir4m e <* 

E c n e e 3 lnA ' {) 

where n e is the thermal electron density, m e is the electron rest mass, e is the electron 
charge, In A is the Coulomb logarithm, eo is the dielectric constant and c is the speed of 
light. Also, in many studies, the runaway electrons were assumed to be ultra-relativistic 
(velocities within 5% to the speed of light) and simplified resonance conditions were used to 
describe the wave-particle interaction. However, the electric field is not always much larger 
than the critical field and the velocity of the electrons is often not that close to the speed 
of light. An example of this is the observations of superthermal electron populations in the 
T-10 tokamak during magnetic reconnection events, when the electric field was transiently 
larger than the critical field during the reconnection (for about 0.1 ms) but then it dropped 
to values near or even below the critical field [jj. Recent work [l(]] has shown that even 
in disruptions the electric field in the core region of the plasma is only slightly above the 
critical electric field, a^l. 

The purpose of this work is to determine what waves could be destabilized by runaway 
beams in a near-critical field. Investigating the lowest relevant limit of the electric field 
when runaway production occurs is a step toward generalizing the analysis of the runaway 
electron driven instabilities to lower electric fields. This way we can gain confidence that 
the analysis of the wave-particle interaction yields valid results in both the high electric field 
and the near-critical limit, before proceeding to the numerical analysis of the interaction for 
electric fields in between. 



In the present work we use the runaway distribution derived in Reference ll| , appropriate 
for a near-critical field, to calculate the instability growth rate of these waves and determine 
the frequencies and wave numbers of the most unstable waves for various parameters. We 
use a general resonance condition, without the ultra-relativistic assumption, so the model 
can be applied also for electrons with lower energies. We show that the whistler branch 
is destabilized via the anomalous Doppler and Cherenkov resonances. Increasing magnetic 
field leads to increasing wave number and frequency while decreasing propagation angle for 
the most unstable wave. The observation of these waves could help to determine the origin 
and evolution of the energetic electrons. If the waves grow to significant amplitude they 
may contribute to efficient transport of particles out from the plasma. 

The remainder of the paper is organized as follows. In Sec. [Til the wave dispersion equation 
is presented, together with a perturbative approximation of the instability growth rate. In 
Sec. IIHI the runaway electron distribution in a near-critical electric field is analyzed and the 
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runaway contribution to the susceptibilities is calculated. In Sec. IIVI the instability growth 
rate of the high-frequency electromagnetic waves driven by runaways is calculated and the 
parameters of the most unstable wave are determined. Here, we also show the stability 
thresholds of the waves and study their parametric dependences. Finally, the results are 
summarized and discussed in Sec. |V] 



II. DISPERSION RELATION 



The dispersion relation of high frequency electromagnetic waves is given by 12 1 



(e n - k\c 2 /uj 2 ) (e 22 - k 2 c 2 /u 2 ) + e\ 2 = 0, (2) 

where u is the wave frequency, k the wave number and e the dielectric tensor of the plasma. 
Equation (T5]) follows from the wave-equation, with the approximation 633 ^> n 2 cos6 l sin^, 
where n = kc/cj is a dimensionless vector with the magnitude of the refractive index, cos 9 = 
k\\/k, 9 is the pitch angle. The subscripts n and j_ denote the parallel and perpendicular 
directions with respect to the magnetic field. The dielectric tensor is 

e = l + X l + X e + X r , (3) 

where x s is the susceptibility of plasma species s, where i denotes the ion, e the thermal 
electron and r the runaway electron populations, 1 is the dyadic unit. As the contribution of 
the runaway population is expected to be small, we consider the dispersion of high frequency 
electromagnetic waves without the runaway term and use the cold plasma approximation 
[12J for the ion and electron populations. The contribution of runaway electrons is added 
later as a perturbation, which is justified in the present case since the runaway electron 
density is much smaller than the thermal electron and ion density. 



Electron-whistler wave 



For the frequency range u ce \jrnjrrii <C u, the background ion and electron contributions 



to e are 



12l 



2 2 

e+ i _ e+l _ U pe _ UJ pe UJce , * 

e ll — e 22 — 1 ^ T aIia 6 12 — 1 I 2V ^ ' 

Here u pe and u ce are the electron plasma and cyclotron frequencies, respectively. Without 
runaways, the dispersion relation can be written as 



£{u) = u\2u? ve + u 2 ce + (k 2 + k 2 )c 2 ] 
WK e + (k 2 + k 2 )c 2 {u 2 pe + u 2 ce ) + k 2 kY\ - k 2 k 2 cW ce = 0. (5) 



Equation (jSJ) has three solutions for u 2 and these can be determined analytically, although 
their closed form expressions are very complicated. One of the solutions satisfies to < knc 
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for all wave numbers k and propagation angles 9 and will be called 'electron-whistler' wave, 
because in certain limits, as we will show, its dispersion characteristics are the same as the 
whistler wave's. For the two other solutions u > k\\c is satisfied. For typical experimental 
parameters, the solution of the analytical dispersion relation (JSJ) has excellent agreement with 
the numerical solution of the full dispersion relation using the hot plasma susceptibilities 
for both ions and electrons from Reference 12] instead of Equation (Tj0). Figure [Th shows 
the three solutions of Equation (]5]) together with the solution of the numerical dispersion 
relation. The solution for the wave frequency is plotted as function of wave number for 
propagation angle 9 = tt/6. The wave dispersions in Figured] are calculated for T = 20 keV. 
The agreement is even better at lower temperatures. 
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FIG. 1: (a) Solution of the analytical approximation of the dispersion relation from Equation ([5]) 
(solid) together with the numerical solution using the hot plasma susceptibilities, for plasma tem- 
perature T = 20 keV, density n e = rtj = 5 • 10 19 m~ 3 , magnetic field B = 2 T and propagation 
angle 9 = ir/6. Dashed line shows to = k\\c. For the electron- whistler wave to < fcyc. (b) The lowest 
frequency solution of the analytical approximation of the dispersion relation Equation §5§ (blue 
solid) together with the whistler approximation from Equation ((6]) (red dashed), for propagation 
angles 8 = (thick lines) and 9 = tt/3 (thin lines). 



The whistler approximation is usually defined by 13] 



fc2c2/ ^cos^-l) = ^|. (6) 



UJ 2 V CO J CO 2 



To show the whistler character of the lowest frequency solution of S(u) = 0, we plot it 
together with the solution of Equation ([6]) as functions of k for 6 = and 9 = tt/3, see 
Figure [Tb. For 9 = there is very good agreement between the two solutions. For 9 = tt/3 
and wave numbers up to 1000 m _1 the two solutions overlap, but for higher wave numbers 
they deviate. This difference is due to the approximation e 33 3> n 2 cos 9 sin 9 used when 
deriving Equatio n flp , while when deriving the whistler wave dispersion given in Equation 
OH]) in Reference [13( no such approximation was used. By investigating the validity of the 



5 



(JSJ) dispersion we concluded that it yields valid results compared to the general dispersion 
relation for magnetic fields up to 3 T. In the following we will therefore limit our analysis 
to B < 3 T. 

Including runaways, Equation (JSJ) can be written as 

where xlj denotes the runaway contribution to the susceptibility tensor. The linear growth 
rate of a small perturbation of the wave frequency u = ojq + 8oj, is 7$ = ^sSco and is given by 





All I w 2 t 22 J 


+ Xh ~ e n) - 2e? 2 Xi2 




2 {3oj$ - 2uj 2 [2uj 2 pe + co 2 ce + (k 2 + kl)c 2 } ~ 


H<+(Ar» + A|)c a « + a4) 





(8) 

where 9 denotes the imaginary part and e°- are the cold plasma dielectric tensor elements 
defined by Equation (TjJ evaluated at the unperturbed wave frequency: e°- = e? +i (u; = ojq). 
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(7) 



B. Magnetosonic-whistler 

To evaluate the wave-particle interaction in the lower frequency region we analyze the 
dispersion relation in the frequency range u ci <w< u ce . Interaction between these waves 
and strongly relativistic runaways has been studied before 0, However, in a near-critical 
field, it is more likely that the runaways are mildly relativistic, and as both the distribution 
function and the resonance condition is different, the analysis in previous work has to be 
generalized. In this frequency range the contributions to the dielectric tensor elements are 

2 ^ 

^ = 1-4 + — ( 9 ) 

CO U c iU ce 

1 1 2 

e+i _ . "pi 

tin — , 
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where u P i and u C i are the ion plasma and cyclotron frequencies, respectively. Substituting 
these into Equation ([2]) leads to the following dispersion relation 



2„,2 



k\v\ k 2 \ ( (k 2 + k 2 -2u 2 /c 2 )v 2 A (k 2 + k 2 )v 2 A tf A tf 




= M (w) = 0, (10) 



6 



where va = cui C i/uj p i is the Alfven speed. In Reference |5|] a simplified version of this 
dispersion relation 



k 2 v\ | 1 4 



(k 2 + k 2 )v\ 



0J c iUJ c 



M s {u) 



0. 



(11) 



valid in the limit u 2 <C k 2 c 2 , was used to study destabilization of waves by an avalanching 
runaway electron distribution. The wave determined by Equation ( fTOj) can be identified 
as the generalized magnetosonic- whistler wave, as its simplified limit, Equation (TIT]) , for 
quasi-perpendicular propagation \k\ 3> \k\\ \ and k 2 c 2 <C ui, 



2J2 



k 2 v\ I 1 + 
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0. 
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(12) 
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has been previously identified as the magnetosonic- whistler wave [5[. 
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FIG. 2: (a) Comparison of the lowest frequency solution of Equation ([5]) (blue solid) with the 
magnetosonic-whistler wave of Equation (|10p (red dashed). The parameters are the same as in 
Figure Q3 (b) Contour plot of the lowest frequency solution of Equation ([5]) and the solution of 
Equation ([6]). The values plotted are u/u: ce on both figures. 



Figure |2] shows contour plots of the electron-whistler wave together with the 
magnetosonic-whistler wave (Figure Eh) and the electron-whistler wave together with 
the whistler approximation from Equation fl6]) (Figure |2b). The electron- whistler and 
magnetosonic-whistler waves approximately overlap with the whistler approximation for 
low k (any propagation angle) and quasi-perpendicular propagation (any k). In the rest of 
the k-6-space the electron-whistler and magnetosonic-whistler waves have different disper- 
sion characteristics, as the magnetosonic-whistler approximation is not valid in the region 
of very high k numbers because its frequency is assumed to be to <C u ce . 
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Including runaways, Equation (ITU]) can be written as 

,2 



M(uj) 
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and the linear growth rate of a perturbation of the wave frequency is 
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"li + ( fc2 + k l - 4cu 2 /c 2 )c 2 (a; ci M e ) + (H + k 2 - 2uj 2 /c 2 )v 2 a 



(13) 

In the following section we will calculate the runaway contribution to the susceptibilities 
which will allow us to evaluate the linear growth rate of the wave. 

III. RUNAWAY CONTRIBUTION 



The susceptibility due to the runaway electron population is given by |X2j ] : 

,2 1*00 poo Cj 

] / 2-Kp ± dpt I dp\\- — 

J J -00 



X 



UJUJn 



ijj 



k\\v\\ 
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where 



m 2 J 2 tt 
—T nL P±U 



. im J^k p±u (j' m y p± u 



u = ?Il + h. ( v ?Il- v ?Il 

dp± uj \ ± dp\\ " dpj_ 



(14) 



Q e = uJce/'j is the relativistic cyclotron frequency of the electrons, J m (z) is the Bessel function 
of the first kind, J' m {z) = dJ m /dz, z = k_ L v±/Q e = k±cp±/uj ce , p = jv/c is the normalized 
relativistic momentum, 7 = y'l + p 2 is the relativistic factor, f r = f/n r is the normalized 
runaway distribution and m is the order of resonance. The general (and implicit) condition 
for the resonant momentum is 

w 7 - muj ce 
n = he ■ 



(15) 



If the distribution function is known, the resonance condition allows the integral in ( 114)) to 
be evaluated using the Landau prescription. 



A. Distribution of the runaway electrons 

To calculate the runaway susceptibilities, the runaway distribution given in Equation (83) 
of Reference [111 is used for the near-critical a ^ 1 case 



tt 1 A ( (° + \ w f, c - , (° + M ^ 

AO-i.pJ = p ( C ,-„/(»-i) °i» {-WTz^y 1 [ l - —i^WTzWt) • (16) 



s 



where 



Z is the effective ion charge and \F\ is the confluent hypergeometric (Kummer) function. 
The distribution function given above was obtained by matching asymptotic expansions in 
five separate regions in momentum space. The calculation is similar to the one presented 



by Connor and Hastie [14J of runaway electron generation, but it is valid for near-critical 
electric field. Note, that to have a positive distribution function, the first argument of \F\ 
should be positive, leading to the condition 1 > C s /(a + 1). Furthermore, the condition 
f r — > as pi: — > oo requires that C s > 2. This gives a region in the a-Z space where 
Equation ( II 6p is valid. The parameter C s as function of a and Z is plotted on Figure El The 
region between the solid and dashed lines gives the combinations of a and Z for which the 
condition 2 < C s < 1 + a is fulfilled. This gives a restriction on the effective charge number, 
since if a ~ 1, the charge number can only be slightly more than unity. In tokamak plasmas 
Z seldom exceeds values of about 3. In the following we will only consider combinations of 
a and Z such that 2 < C s < 1 + a. One such combination is a = 1.3 and Z = 1 and this, 
together with the parameters n e — 5 • 10 19 m~ 3 , B — 2 T, are the baseline parameters of 
our study, and will be used in the rest of the paper unless otherwise is stated. Note, that 
although C s includes a term proportional to \j\Ja — 1, its value varies very little in the 
parameter space where the distribution function is valid, as Figure [3] shows C s is between 
2.5 and 3 in the region of interest, irrespective of the exact value of a and Z. Therefore the 
distribution function is not very sensitive to these values. 

Equation (Tl6l) is valid for all p > p c in the case of near-critical electric field. Note, that 
the integral of Equation (fT6l) function in the whole momentum space is divergent. This is 
because the electric field continuously accelerates electrons and more and more electrons 
will run away. In spite of the continuous acceleration, the distribution is in quasi-steady 



state, as the water leaking out of an unplugged bath tub [15j. However, as the existence 
of the electric field is finite in time, there is a maximum number of runaways and there 
is a maximum energy which runaway electrons can reach in reality. In the expressions for 
the runaway susceptibilities we use a normalized distribution function J f r d 3 p = 1. The 
normalization constant A in Equation ( 1161) is obtained from 



Pmax 



dp±2np ± dp\\f r (p\\,p±) = l, (18) 

Jp c 

where p max is the normalized momentum corresponding to the maximum energy. This 
integral can be easily solved numerically if p max is known. The value of p max depends on 
the exact value and time evolution of the accelerating field. In this paper we approximated 
the maximum energy as 2.6 MeV, corresponding to J5 max = 5. A typical value of the 
perpendicular momentum can be determined from the runaway distribution function. In 
the case of p\\ m ax = 5, this is p± ~ 3. This value corresponds to E± = 1.6 MeV. Figure 0] 



9 




1 2 3 4 5 6 7 8 
Z 



FIG. 3: C s as function of a and Z . The distribution function is valid in the region 2 < C s < 1 + a. 
Solid black line shows C s = 1 + a and dashed black line is C s = 2. The region between the solid 
and dashed lines gives the combinations of a and Z for which the condition 2 < C s < 1 + a is 
fulfilled. 

shows Equation (TTB]) for Z — 1 and a = 1.3. For the baseline parameters of our study this 
corresponds to the electric field of 0.06 V/m, while the critical field is 0.046 V/m. 



a = 1.3 




20 2 

FIG. 4: Normalized runaway electron distribution function in near-critical field, f r /A plotted with 
respect to the parallel and perpendicular momentum normalized to m e c, for Z = 1 and a = 1.3. 



It is instructive to compare the distribution in Equation ( TT61) with the distribution derived 
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for the case of secondary runaway generation 

where a = (a — 1)/(Z + 1) and cz = \/3(Z + 5)/7rln A. The avalanche distribution is based 
on the solution of the kinetic equation for relativistic electrons in the limit of a ^> 1 using 



the Rosenbluth-Putvinski runaway growth rate 161 ] 

dn r n r (a — 



dt CzT 

as boundary condition. Here r is the collision time for relativistic electrons. This means 
that the runaway density grows exponentially as n r = n r oexp [(a — l)t/(rcz)], where n r0 
is the seed produced by primary generation. Note that for the avalanching distribution 
2tt J p±dp±dp\\f r = 1, independent of the maximum momentum. The distribution (fT9|) is 
valid if secondary generation of runaways is dominant, as expected to be the case in large 
tokamak disruptions. In contrast, the distribution 016p is valid when primary runaway 
production is the main source of the superthermal electron population. 

Comparing the two distributions, we note that the near-critical distribution function in 
Equation ( jl6]l represents a broader beam, with a less rapidly decaying tail. Figures [ME] show 
the comparison between the near-critical and the avalanching distribution functions. Figure 
[5^ shows the near-critical distribution for two different values of a and Z = 1.5. Figure 
[5b shows the comparison between the a = 1.3, Z = 1 case of the near-critical distribution 
with the avalanching distribution, which is significantly more beamlike. To illustrate that 
the distribution is more beamlike and more rapidly decaying in p\\ in the avalanche case, in 
Figure |6] we show the comparison between the near-critical and avalanching distributions for 
specific values of p\\ and p±. In spite of the differences noted above, the distribution functions 
in the a £ 1 and a ^> 1 limits are similar in the sense that both have an anisotropy in the 
p\\ direction, and a smooth transition between the two can be envisioned based on Figure |5j 
The reason for using this particular distribution function (Eq. (fT6|) ) in the present work is 
that it is at the lower limit of a that can possibly produce runaway electrons. 

Although Equation ffTB]) describes primary generation of runaways, it does not mean that 
the generation rate is small. Primary generation implies that the runaway generation is 
smaller than n e /r. But since n e /r is very large, primary generation can result in a sub- 
stantial runaway electron population and its importance has been shown in many numerical 



simulations, see e.g. [H| 



B. Resonance condition 



In a plasma with a slightly supercritical electric field, the characteristic value of the 
normalized momentum p in the runaway region satisfies p > 1/ y/a — 1. To obtain an 
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FIG. 5: (a) Contour plot of the distribution function, f r /A for a = 1.3 (solid, corresponding 
to E = 0.06 V/m) and a = 1.5 (dashed, E = 0.069 V/m). The effective charge is Z = 1.5. 
(b) Comparison between the near-critical, f r /A (blue solid) and avalanche, 100/ r dlsr (red dashed) 
distribution functions. For the near-critical distribution we used Z = 1 and a = 1.3. For the 
avalanche distribution we used InA = 18, Z = 1 and E 1 = 40 V/m (corresponding to a = 865). 



explicit formula for the resonant momentum, the expression 7 = \/^- + P 2 ± + P 2 should be 
substituted into the resonance condition, and that leads to 



—k\\ cmu ce ± UqJ (kjc 2 — Wq)(1 + p 2 ± ) + m 2 u 2 e 

Pros (p_L, fc|| , Wo) = fc 2 c 2 _ ^2 " ( 20 ) 

By using this general resonance condition, the expressions giving the imaginary parts of the 
runaway susceptibilities become quite complicated. The full expressions for the susceptibil- 
ities are given in Appendix A. 

Only the p rcs > resonant momenta are physically relevant. By studying the p res > 
condition for different signs of m, using the relation between k\\c and uo{k,8) it can be 
shown that the Doppler resonances (m > 0) cannot be satisfied for any of the solutions in 
Equation (jSJ) or Equation ( fTOl) . 

a. Anomalous Doppler resonance For the anomalous Doppler resonance (m < 0) the 
p res > condition, defining the physically relevant region of the p res resonant momentum, is 



k\\c\m\uj ce + uo x j (k 2 c 2 - Wq)(1 +p±) + m 2 u 2 e 

^ > 0, (21) 



k 2 c 2 



leading to k 2 c 2 > ^(/c, 0), which is only satisfied for the electron- whistler branch and not 
for the other two solutions of Equation fl5]). Also the magnetosonic- whistler wave can be 
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FIG. 6: Comparison between the near-critical, f r /A (blue solid) and avalanche, 100/^ lsr (red 
dashed) distribution functions. For the near-critical distribution we used Z = 1 and a = 1.3. For 
the avalanche distribution we used In A = 18, Z = 1 and E = 40 V/m (corresponding to a = 865). 

(a) The distribution function as a function of pu for p± = (thin lines) and p± = 0.5 (thick lines). 

(b) The distribution function as a function of p± for pu = 20 (thick lines) and jpn = 40 (thin lines). 



destabilized via this resonance. 

b. Cherenkov resonance For Cherenkov resonance (the case of m = 0) the p vcs > 
condition is 



uj 0x l{ky-ul){l+pl) 

> 0. (22) 



He 2 



ijj, 



2 



^0 

This also leads to the condition k^c 2 > u 2 (k,9), narrowing down the possible waves once 
again to the electron-whistler wave and the magnetosonic-whistler wave. Summarizing the 
results above, we conclude that for m < 0, only the electron-whistler waves can yield phys- 
ically relevant results out of the high frequency electron waves denned by the dispersion 
relation in Equation (J5J). The magnetosonic-whistler waves can also be destabilized via 
m < resonances. However, combining the region of the validity of the wave frequency 
with the resonance condition it can be seen that in the magnetosonic-whistler case the 
destabilization is most effective by very energetic (around 10 MeV) runaway electrons. 

If p ^> 1, for the beam-like distribution function in Equation (Tl6l) with p\\ ^> p±, the 
7 ~ \p\\\ approximation can be used (which will be called the ultra-relativistic limit), and 
the resonance condition in (fl4"]) simplifies to 

-muj ce 

P\\ = 1 23 

k\\c — oj 



and m < for physically relevant results. 
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FIG. 7: Normalized growth rate 10 3 7i/w ce for the electron- whistler wave (a,b) and the 
magnetosonic- whistler wave (c,d). Both in (a,b) and (c,d) the black line is co = w ce /45, the 
electron-whistler approximation is valid in the region above it. In (c,d) the dashed line denotes 
(jj = uj c i, the magnetosonic- whistler approximation is valid in the region above it. The rest of the 
parameters are n e = 5 • 10 19 m -3 , n r = 3 • 10 17 m~ 3 , B = 2 T and p mSLX = 5. (a,c) Ultrarelativistic 
resonance condition for m = —1. (b,d) General resonance condition, sum of the cases m = —1 and 
m = 0. 



IV. UNSTABLE WAVES 

The instability growth rates for the electron-whistler and the magnetosonic-whistler waves 
can be calculated from Eqs. OH]) and ffl3|) as functions of k. Figure [T^b shows the growth 
rates for the electron-whistler wave using the ultrarelativistic limit and the general reso- 
nance condition, respectively, for the baseline parameters. The growth rate increases with 
decreasing k throughout the range of validity of the electron-whistler approximation. Com- 
paring Figure [T^l and Figure [TJd it can be seen that by using the ultrarelativistic condition 
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one gets somewhat different results than by using the general condition, but the qualitative 
behaviour is the same. 

Figure [Tbd shows the growth rates for magnetosonic-whistler wave. In contrast to the 
electron-whistler wave, the growth rate has a well-defined maximum. However, as it was 
mentioned before, the resonance condition cannot be satisfied in the region close to the maxi- 
mum growth rate unless the resonant energy is very high (around 10 MeV for the parameters 
given in the figure caption), which is only expected to be reached by few electrons in the 
near-critical case considered in this paper. Interestingly in both cases, the largest instability 
growth rate occurs for the region in the fc-#-space where the whistler approximation fl6]) is 
valid (see the low k and perpendicular propagation part of Figure [2]). However, we note that 
only high energy electrons can interact with the quasi-perpendicularly propagating whistler 
wave, therefore it is more likely that the electron-whistler branch with slightly higher k and 
more oblique propagation is the one which is the most unstable wave. 

A. Most unstable wave 

The normalized momentum corresponding to the maximum energy of 2.6 MeV is ap- 
proximately p TCS = 5. The corresponding wave numbers and propagation angles of the 
electron-whistler wave can be calculated by using the general resonance condition and the 
dispersion relation. The growth rate and the values corresponding to p rcs = 5 are shown on 
Figure El The most unstable wave in the near-critical case is an electron-whistler wave with 
frequency 4.2- 10 10 s" 1 ~ 0.12u ce (for a magnetic field of 2 T), wave number of approximately 
650 m _1 , and angle of propagation 9 « 0.9. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
6 [rad] 



FIG. 8: Most unstable wave in the near-critical case: maximum of the growth rate (10 3 7i/a; ce , 
contour lines) on the line corresponding to the maximum runaway energy (2.6 MeV, white dots). 
The parameters are n e = 5 • 10 19 m -3 , n r = 3 • 10 17 m~ 3 , B = 2 T. 
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It should be noted, that the parameters of the most unstable wave are sensitive to the 
magnetic field. The reason is that the resonance condition is highly dependent on the 
magnetic field through the gyrofrequency. Due to this fact, the p < 5 condition for the 
momentum of the runaway electrons yields very different wave numbers for the most unstable 
wave. For example, if the magnetic field is 4 T instead of the 2 T on Figure [U the p Tes = 5 
resonant momentum yields k ~ 1600 m _1 wave number and 9 ~ 0.3. Therefore, by increasing 
the magnetic field, the wave number and frequency of the most unstable wave increase, while 
the angle of propagation decreases. 

The wave number, propagation angle and frequency of the most unstable wave also depend 
on the maximum runaway energy. Figure [9] shows that as the energy grows the propagation 
angle becomes larger and the wave number and frequency drops. This means that for low 
energy runaway electrons (energies just above the critical energy for runaway acceleration) 
we expect frequencies slightly less than half of the electron cyclotron frequency, propagating 
angles of 9 ~ 0.4 and wave numbers of 1000 m _1 . As the runaway energy grows the frequency 
and wave number of the most unstable wave falls. 
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FIG. 9: The value of wave number (blue dashed) and propagation angle (red dotted) (a) and 
frequency (b) of the most unstable wave as function of maximum runaway energy. 
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B. Stability diagram 

In order to determine the stability limits, the instability growth rate of the wave has to be 
compared to the damping rates. In cold plasmas, collisional damping is dominant, and the 



damping is approximately equal to = l-5 T ei 18(] , where r ei = 37r 3 / 2 mg t4 e eg/7ijZ 2 e 4 In A 
is the electron-ion collision time. In addition to collisional damping, the wave is damped due 
to the fact that the extent of the runaway beam is finite, and the wave energy is transported 
out of its region with a du/dkj_ perpendicular group velocity. This mechanism can be 
accounted for by adding a convective damping term j v = (doj / dk±) / (AL r ) , where L r is the 



radius of the runaway beam [6j. The linear growth rate of a wave is thus 7; = 7« — 7d — 7^, 
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and the wave is unstable, if 7^ > 0. A simple estimate of the order of magnitude of the 
damping rates shows that for typical parameters, and for reasonably narrow electron beams, 
the convective damping is expected to dominate if T e > 200 eV. However, as the plasma 
temperature seldom reaches such a high value in the relevant case of tokamak disruptions 



19[, collisional damping should not be neglected. 




1.5 2.0 

B[T] 



FIG. 10: Stability thresholds for the most unstable wave in near-critical electric field, for electron 
temperature T e = 20 eV. (a,b) Stability threshold as function of magnetic field for the electron- 
whistler wave and magnetosonic- whistler waves, respectively. The runaway-beam radius is L r = 0.1 
m (dashed) and L r = 0.2 m (solid). In (a) we assume p max = 5 and in (b) p max = 20. (c,d) 
Sensitivity of the stability threshold to the normalized electric field a for the electron-whistler 
wave. The runaway beam radius is L r = 0.1 m, n e = 5 • 10 19 m~ 3 and the maximum runaway 
energy is 2.6 MeV, corresponding to p m ax = 

5. In (c) Z = 1 and in (d) Z = 1.5. 



The linear stability threshold was determined in the following way. For any given value 
of the magnetic field the growth rate is calculated (for m — — 1, m — 0, then adding them 
for all k and 8), then the collisional and convective damping rates are subtracted from it. 
The parameters of the most unstable wave are then determined. The stability threshold 
of the most unstable electron-whistler wave is shown on Figure [T0a . We conclude that for 
typical parameters the runaway density needed to counter the damping rates, therefore to 
destabilize an electron-whistler wave is of the order of 10 m -3 or n r /n e = 0.2%. For the 
magnetosonic-whistler wave, the stability threshold is shown in Figure [TOb . In this case we 
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assumed p max = 20, since in this case the resonant particles have higher energies than in the 
electron-whistler case. Figure [TUbd shows the dependence of the stability threshold on the 
normalized electric field a. We conclude that above B >, 1.5, the runaway density needed 
for destabilization is sensitive to a, for higher normalized electric field lower runaway density 
is needed to destabilize the wave. This dependence on a is due to the fact that for a higher 
electric field the anisotropy of the runaway distribution and thus the destabilizing effect is 
stronger, therefore a lower density of runaways suffices for a resonant destabilization of the 
wave. This is by no means because of the special characteristics of the model distribution 
we used, but is due to the underlying physics. 

It is instructive to compare the order of magnitude of the runaway density required for 
destabilization of the whistler wave to the one that is measured in an experimental setup. 
Figure [11] shows the stability thresholds as function of magnetic field, for various runaway 
beam radii, for parameters relevant to an experiment in the T-10 tokamak j^]. We note that 
the runaway density needed for destabilization is about 10 17 m -3 even for a narrow runaway 
beam. The runaway density estimated in the experiment was almost an order of magnitude 
higher than this: 7 • 10 17 m~ 3 . 
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FIG. 11: Stability threshold for the most unstable electron- whistler wave in near-critical electric 
field, for the experimental parameters of the T-10 tokamak. The parameters are a = 1.9, Z = 3, 
n e = 4 • 10 19 m- 3 , T e = 0.5 keV, p max = 1.5. 




V. CONCLUSIONS 

The presence of high energy electrons is often associated with bursts of high-frequency 
waves. The emission of radiation is most often due to Bremsstrahlung and synchrotron 
radiation, but in certain cases they are due to instabilities caused by the anisotropy. The 
observation of these waves can help to determine the origin and evolution of the energetic 



electrons, and also in some cases, the properties of the background plasma 13] . The insta- 



bility may result in pitch-angle scattering induced isotropization and may therefore prevent 
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the harmful effects of the runaway electron beam 

The reason for the generation of an anisotropic runaway electron population is the high 
electric field that is often caused by reconnection events in magnetized plasmas. In previous 
calculations regarding waves driven by runaways the electric field was assumed to be much 
higher than the critical field a>l. This is not often the case in reality. Therefore, in this 
paper we use an electron distribution function that is valid in the near-critical case. We 
show that in this case, the distribution is broader and less rapidly decaying compared to the 
a ^> 1 case. 

By studying the linear growth rate of the electron- whistler branch (valid in the frequency 
region Lo ce \fm^rni <C u) and the magnetosonic- whistler branch (valid in the frequency 
region UaCuC u ce ) separately, we find that the frequency of the most unstable wave is in 
the region where these overlap and have characteristics similar to the whistler approximation. 
For typical tokamak parameters we find that the frequency of the most unstable wave is 
around 0.12u ce , in the near-critical case with a = 1.3 and E = 2.6 MeV. The frequency 
and wave number of the most unstable wave depends strongly on the magnetic field and on 
the maximum runaway energy. By comparing the ultra-relativistic limit of the resonance 
condition and the general one we show that although the behaviour of the instability growth 
rates of the electron-whistler and magnetosonic-whistler waves are similar, the actual values 
for the growth rate may differ, and therefore the frequency and wave number of the most 
unstable wave might be different. 

The instability growth rate of the electron-whistler wave was compared to the collisional 
and convective damping rates. We find that the number density of runaways that is required 
to destabilize the waves increases with increasing magnetic field. For low magnetic fields the 
convective damping decreases, while the collisional damping rate remains constant, making 
it dominant in this region. As the growth rate also decreases, the stability limit is high 
for low magnetic fields. We investigated the stability of the whistler waves for parameters 
relevant to the T-10 tokamak {{J, where the effective electric field is near-critical. We found 
that the observed runaway density is about an order of magnitude higher than the density 
needed for the most unstable electron-whistler wave to be destabilized. Thus the runaway 
population may indeed give rise to this whistler wave. 

The importance of this study is that it considers the case where the electric field is near 
critical, which is opposite to the other limit that has been considered in previous work 0- 
3 (when the electric field is far above the critical). By investigating this case, we show 
that the high-frequency instabilities are qualitatively similar, but have different frequencies 
and wave numbers. This result may open up the possibility of diagnostics. Understanding 
the properties of the waves destabilized by runaway electrons can be important in view of 
obtaining information about the energetic electron population and the background plasma. 
Regardless of the fact that the distribution used in this work is only valid in the near- 
critical case, if we compare it to the avalanche distribution we observe a smooth transition 
between the two, and so we expect that the distribution does not change qualitatively. Also, 
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the characteristics of the growth rates in the near-critical and high electric field limits are 
similar, in the sense that the maximum of the growth rate is at low wave numbers and 
near-perpendicular propagation in both cases & The d.fferences in the parameters of the 
most unstable wave for near-critical and avalanching cases are mainly due to the maximum 
runaway energy. The similarity of the results, added to the relaxation of the approximations 
used in previous work, opens the way toward more general numerical studies of wave-particle 
interaction for arbitrary electric fields. 

The whistler waves, if they grow to significant amplitude, in principle could perturb the 
background magnetic field and lead to efficient transport of particles (specially energetic 
ones) out from the plasma. The effect of magnetic field perturbation has been studied 



before 2CH22j . and it has been shown that runaway avalanches can be prevented altogether 



with sufficiently strong radial diffusion. However, the magnetic fluctuation level that is 
required for this to happen is estimated to be 5B/B ~ 10 -3 and the magnetic fluctuation 
level induced by these high-frequency whistler waves would be several orders of magnitude 
lower than this value. Therefore, as mentioned before, the runaway electron population 
will be affected mostly through pitch-angle scattering and concomitant isotropization and 
synchrotron radiation damping. 
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Appendix A: Runaway susceptibilities 

In the general case the susceptibilities have the following form: 

Imx^k,^) = nflT / d P± / d Pll Y,rn 2 J 2 m (z) x (24) 

( df r (p) ( mu ce \ df r (p) k\\cp± \ 1 ^ ( _ fcjjggjj _ mu ce \ 
V dp± V 7 / 9p\\ 7 / 7 V ° 7 7 / 

where we used the resonance condition 7^0 — kucp\\ = mu ce to replace the factor 
(u — k\\cp\\/'~i) . The Inr^ an d IniXi2 terms only differ in multiplicative constants and 
will be presented later. 



20 



For a general function g(p\\) we can rewrite the integral in p\\ as follows 



7 



1 



5x D 
dx5 ( A , + 



<?(x). (25) 



where a; = p\\, A — ujq, B = k\\c, C = 1 + p 2 ± and D = —mu) ce . Changing variables 
y = (Bx - D)/y/C + x 2 , so that 

dX BC + xD V ' 

BD ± yy/(B 2 - y 2 )C + D 2 



x 



B 2 -y 2 



Equation (l25j) yields 

(C + x 2 ) 3 / 2 



dy 



8{A-y)g 



BD±yy/(B 2 -y 2 )C + D 2 



BC + xD v ' \ B 2 -y 2 
(AD ± B^(B 2 - A 2 )C + D 2 ^ ' 



(26) 



( BD± Ay/(B 2 — A 2 )C + D 2 



(B 2 - A 2 ) 3 



BC 



D 



B 2 —A 2 



BD ± Ay/ (B 2 — A 2 )C + D 2 ) \ 



B 2 -A 2 



Using the above expression to solve the integrals in the runaway susceptibilities, we arrive 
to the following formulas: 



ImXn (k, UJo) 



27i 2 uj 2 pr u 2 ce 



ulk\c 2 



poo 

/ (ij)iVm 2 ^(z)x 
Jo 



(27) 



df r (p) (mu ce \ df r (-p)k\\cp L 



dpi 



7 



dp\\ 7 



h (p±,k\\,u ) 



P\\=Pr. 
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where 

h (p±,k\\,u ) 



1 



(k 2 c 2 - w 2 ) 3 

-uj mujce ± fc||cy (^ 2 c 2 - w o) (1 + P±) + m2uJ cej 



(28) 



fc||c(l+^i) - (^-A;|| C ma; ce ± (*fc 2 - ^ 2 ) (1 + ^) + m 2 u 2 c 



Similarly the other two terms of the runaway susceptibility: 



ImX22 ( k > w o) 



2n 2 u 2 pr 



ijj, 



o Jo 



P idp ± j2(j^)y 



X 



9/r(p) (mu ce \ 8 f r {p) k\\Cp± 



dpi 



7 



9p|| 7 



7 



(29) 
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Rex^ 2 (k, wo) 



27r 2 u;;Lc<; r 



3/r(p) 

dp± 



poo 

/ Px^P-l ^ mJ,,. , 
Jo 



2)^(2) x 



(30) 



7 



dp,. 







/i (p_L, fc||,o;o) 



P||=Pr 



7 



Equation f )30|) gives the real part of the runaway susceptibility X12 since this term is the one 
needed in the expression for the growth rate. By substituting these runaway susceptibilities 
into the expression of the growth rate, the calculations yield results in the general, relativistic 
case regarding the runaway electrons interacting with the corresponding wave. 
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